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ABSTRACT 

We identify two candidate magnetars in archival X-ray observations of HESS detected shell-type 
SNRs. X-ray point sources in CTB 37B coincident with HESS J1713-381 and in G353.6-0.7 coinci- 
dent with HESS J 173 1—347 both have AXP-like spectra, much softer than those of ordinary, rotation 
powered pulsars, and no optical/PR counterparts. The spectrum of CXOU J171405. 7— 381031 in 
CTB 37B has a hard excess above 6 keV, which may be similar to such components seen in some 
AXPs. A new Chandra observation of this object reveals a highly significant pulsed signal at P = 3.82 s 
with pulsed fraction f p = 0.31. Analysis of an XMM-Newton observation of the second candidate, 
XMMU J173203.3-344518 in G353. 6-0.7, yields only marginal evidence for a 1 s period. If it is not 
a magnetar, then it could be a weakly magnetized central compact object (CCO). Considering that 
these HESS sources previously attributed to the SNR shells are possibly centrally peaked, we hy- 
pothesize that their pulsars may contribute to diffuse TeV emission. These identifications potentially 
double the number of magnetar/SNR associations in the Galaxy, and can be used to investigate the 
energetics and asymmetries of the supernovae that give rise to magnetars. 

Subject headings: ISM: individual (CTB 37B, G353.6— 0.7) pulsars: individual 

(CXOU J171405. 7-381031, XMMU J173203. 3-344518) — stars: neutron 



1. INTRODUCTION 

Surveys by the atmospheric Cherenkov telescopes 
HEGRA, HESS, MAGIC, and VERITAS, have re- 
vealed more than 50 Galactic TeV sources 1 , the ma- 
jority of which are supernova remnants (SNRs) or pul- 
sar wind nebulae (PWNe), the latter b eing the largest 
class. For reviews of TeV PW Ne. see iCarrigan et aLl 
(2008 ). IGallant et al.l (l200l. ILemiere et alJ 32008). 
iHessels et al.l (|2008[) . and iMattana et alj (|2009) . In al- 
most all cases, the pulsars responsible for the PWNe 
that are detected at TeV energies have spin-down lumi- 
nosities > 10 36 erg s _1 . Exceptions are the very nearby 
Geminga pulsar (E — 3.2 x 10 34 er g s -1 ), with diffus e 



TeV emission detected by Milagro (jAbdo et ail [ 

and PSR J2032+4127 (E = 2.7 x 10 35 erg s" 1 ) associ- 
ated wit h TeV J2032+4130 in the Cyg OB2 stellar asso- 
ciation (jCamilo et al J 12009) . There is no evidence that 
any of the 14 previously known magnetars 2 produce TeV 
emission. Their spin-down luminosities, E = 4ir 2 IP/P 3 , 
are < 10 35 erg s _1 , less than that of known pulsars pro- 
ducing X-ray PWNe, and it is not even certain that they 
accelerate particles to TeV energies as do ordinary radio 
a nd spin-powered p uls ars. 

lAharonian et al.1 ([2006) discovered the TeV 
source HESS J1713— 381 and associated it with 
the SNR CTB 37B. U sing a Chandra observation, 
lAharonian et al.l (|2008af ) then found the point X-ray 
source CXOU J171405. 7-381031 in CTB 37B, and 
considered that it could be a pulsar, albeit with an 
unusu ally soft, non-thermal spectrum. iNakamua et alJ 
(2009) analyzed Chandra and Suzaku spectra of 
CXOU J171405. 7-381031, noting also possible vari- 

1 VHE T-rav Sky Map and Source Catalog, 
I htt p : / / www . mppmu. mpg . de/ ~ Irwagner /sources / index.html 

2~ jjQR/AXF Online Catalog, 

http:/ /www. physics. mcgill.ca/~ pulsar/magnetar/main.html 



ability in flux, which they took to be good evidence 
that it is an anomalous X-ray pulsar (AXP). In 
Section 2 of this paper, we report the discovery of 
pulsations from CXOU J171405. 7-381031 that veri- 
fies this conjecture. We also analyze archival X-ray 
data on another compact source in a HESS detected 
SNR, G353.6-0 7=HES S J1731-347 (|Tian et all [20081 : 
lAharonian et al.l l2008b). showing that it too is a mag- 
netar candidate based on its spectrum and a possible 
periodicity. Considering these new associations, we 
discuss in Section 3 the possibility that an early stage 
of a magnetar may produce a Te V PWN, a preceden t 
being PSR J1846-025 in Kes 75 (|Diannati-AtaJl25ol) . 
We conclude in Section 4 that it is of great interest to 
measure the spin-down rates of these new candidate 
magnetars to test this hypothesis. 

2. X-RAY OBSERVATIONS OF NEW CANDIDATE 
MAGNETARS 

2.1. HESS J1713-381/CTB 37B 

At radi o wavelengths , CT B 37B is a shell « 10' in 
diameter (iKassim et al.lll99ll ). Its distance, estimated 
from H I absorption a nd the Galactic rotation curve, is 
» 8 kpc reen ||200H ). Its age may be only ~ 1500 yr, 
which m ay associate it with the historical supernova 
373 AD (jClark fc Ste phenson| [l975Tl. Following t he dis - 
covery of HESS J1713-381 bv lAharonian et al l 12006). 
further HESS observations by lAharonia n et al.l |2008a) 
localized it to the center of the radio shell of CTB 37B, 
and found that the extent of the TeV source is compatible 
with either a centrally peaked Gaussian of a = 2.'6 ± 0/8 
or a shell of radius 4'-6', consistent with the size of the 
radio remnant. The energy flux of HESS J1713— 381 from 

, corresponding 
where da is the 
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0.2-10 TeV is « 4.2 x KT" erg on- 
to luminosity L 1 = 3.2 x 10 34 c?g 
distance in units of 8 kpc. 
The Chandra ACIS-I observation listed in Table [T] 
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Fig. 1. — HESS TeV sources denoted by dashed ellipses representing their estimated intrinsic extent (Aharonian et al. 2008a, 2008b). 
X-ray point sources (crosses) are magnetar candidates. Left: HESS J1713— 381 with the shell-type radio SNR CTB 37B (contours); MOST 
843 MHz) and X-ray emission (greyscale; Chandra ACIS-I). Right: HESS J1731-347 associated with the shell-type radio SNR G353.6-0.7 
(contours; ATCA 843 MHz) and its X-ray emission (greyscale; XMM-Newton EPIC MOS). The bright, compact radio source SSW of the 
pulsar is extragalactic, based on H I absorption velocities (Tian et al. 2008). The radio source to the NW is an H II region at 3.2 kpc, 
which is presumed to be associated with the SNR. 



TABLE 1 

Log of X-ray Observations of HESS Sources 







Instrument/ 


ObsID/ 


Date 


Exp. 




Source 


Mission 


Mode 


Seq. No. 


(UT) 


00 


Refs. 


J1713-381 


Suzaku 


XIS/Normal 


501007010 


2006 Aug 27 


82,815 


1 




Chandra 


ACIS-I/TE/VF 


6692 


2007 Feb 2 


25,155 


1,2,3 




RXTE 


PCA 


93011-03-01 


2008 Jul 28-31 


89,326 


3 




Chandra 


ACIS-S/CC 


10113 


2009 Jan 25 


30,146 


3 


J1731-347 


Suzaku 


XIS/Normal 


401099010 


2007 Feb 23 


40,619 


3,4 




XMM-Newton 


pn+MOS/FF 


0405680201 


2007 Mar 21 


25,408 


3,4,5 




Chandra 


ACIS-I/TE/VF 


9139 


2008 Apr 28 


29,610 


3,5 




Swift 


XRT/PC 


00037740001 


2009 Feb 4 


1,387 


3 




Swift 


XRT/PC 


00037740002 


2009 Mar 15 


1,432 


3 



References. — (1) Nakamura et al. 2008; (2) Aharonian et al. 2008a; (3) This work; (4) Tian 
et al. 2009; (5) Acero et al. 2009. 



revealed a bright X-ray point sour c e and faint, dif- 
fuse emission that lAharonian et al.l ([2008a) fitted to 
a non-equilibrium ionization model. In Figure [1] we 
show the Chandra image with superposed radio contours 
from the Molonglo Galactic Plane Survey at 843 MHz 
([Green et al.lll99 9), and the la extent of the TeV emis- 
sion. Although only the eastern edge of the radio 
SNR is seen her e, at lower frequen cy a nearly circular 
shell is evident ([Kassim et al.l 119911 ). The point source 
CXOU J171405. 7-381031 lacks an optical or IR coun- 
terpart in the Digitized Sky Survey or in 2MASS, and 
is likely associated with the SNR ba s ed on its location, 
flux, and spectrum. lAharonian et al.l (|2008af ) noted that 
an absorbed power-law fit to the point source, with pho- 
ton index T = 3.3^^ , yields an Nh that is consistent 
with the fit to the diffuse emission. Alternatively, they 
obtained a good fit with a two-temperature blackbody of 
kTx = 0.52 keV and kT 2 = 1.6 keV. There is no strong 
evidence for a PWN associated with the point source. 



Using the Suzaku observation listed in Table [TJ 
INakamua et alj (|2009f ) found nonthermal X-ray emission 
in the southern region of CTB 37B. They also reana- 
lyzed the Chandra spectrum of the point source, find- 
ing a 2-10 keV luminosity of 1.4 x 10 34 d\ ergs s _1 . 
Comparing that with the Suzaku luminosity, which ap- 
peared to be a factor of 1.8 higher, they conjectured that 
variability plus the steep X-ray spectrum indicated that 
CXOU J171405. 7-381031 is an AXP. 

We obtained a new Chandra observation of 
CXOU J171405. 7-381031 on 2009 January 25 us- 
ing ACIS-S in continuous clocking (CC) mode for 30 ks 
to search for the expected pulsations. After barycentric 
correction, we searched for a periodic signal to the 
Nyquist limit of 5.7 ms. We find a highly significant 
signal at a period of 3.82305 ± 0.00002 s, having a single 
broad pulse with pulsed fraction f p = 0.31, defined 
as the fraction of counts above the minimum in the 
light curve. Figure [2] shows the folded light curve in 
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Fig. 2.— Z\ power spectrum of CXOU J171405.7-381031 from 
the Chandra CC-mode observation, and folded light curve (inset) 
of the new 3.82 s pulsar. Normalized count rate is plotted. 



TABLE 2 

Spectral Fits to Chandra ObsID 6692: 
CXOU J171405. 7-381031 



Parameter PL PL + BB 



JV H (10 22 cm- 2 ) 


4 03+0.51 


2.59 ±0.43 


r 


3.10 ±0.32 


-2 0+ 1 - 1 

z ' u -1.0 


kT (keV) 




64+ - 07 


R (km) 




0.82 d 8 


F x {2 - 10 keV) a 


1.8 x 10 -12 


1.8 x 10 -12 


L x {2 - 10 keV) b 


1.4 x 10 34 d| 


1.7 x 10 34 d| 


x 2 » 


1.28(96) 


1.06(94) 



a Absorbed flux in units of erg cm 2 s 1 . 
b Unabsorbed luminosity in units of erg s — 1 . 



the 1-5 keV band, which contains most of the counts. 
The peak power measured using the Z\ statistic is 194, 
corresponding to a vanishing chance probability. There 
is no evidence for energy dependence of the pulse profile. 
The long period is consistent with expectation based on 
the AXP-like spectrum, and falls in the 2-12 s range 
seen for magnetars. It is possible that the true period is 
7.6 s, with two peaks per rotation, but there is no strong 
evidence to support this in the form of distinguishable 
peaks in the current data. 

Listed in Table [I] is a set of archival RXTE ob- 
servations that were pointed 14' from the position of 
CXOU J171405. 7-381031, well within the response of 
the Proportional Counter Array. We searched these data 
for the 3.82 s period, but no significant signal was found. 
Given the Chandra measured spectrum and flux of the 
source (discussed next), the lack of detection by RXTE 
is not surprising. 

We reanalyzed the spectrum from the ACIS-I obser- 
vation, fitting a power-law or blackbody model as be- 
fore. The results are listed in Table O For the power- 
law model, we obtain a photon index T — 3.1, con- 
sistent with the published analyses. However, this fit 
is not acceptable (see Figure [3j, as there is clear ex- 
cess of flux above 6 keV . T his deviation is not see n in 
lAharonian et alj (|2008aD and iNakamua et all (|2009fl , ev- 
idently because they excluded channels above 7 keV. A 
single blackbody is even a worse fit, and a two blackbody 
model requires an unphysical high temperature. Al- 
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Fig. 3.— (a) Spectrum of CXOU J171405.7-381031 from t e 
Chandra ACIS-I/TE observation fitted to a power-law model. Pos- 
itive residuals are evident above 6 keV. (6) The same spectrum fit- 
ted to a power-law plus blackbody model (see Section 2.1). Fitted 
parameters are listed in Table [2] 

though we can get a adequate fit using a power-law plus 
blackbody model, the resulting negative spectral index, 
T = — 2.0^'J is highly unusual, and possibly not mean- 
ingful. This hard X-ray excess may be an extreme case of 
the components that dominate at energie s above 10 keV 
in several AXPs dMolkoy et all l2004t iKuiper et all 
12004 IRevnivtsev et al.ll2004tlden Hartog et al.ll2006D . al- 
though the latter have spectral indices 0.9 < T < 
1.4. A more unusual possibility would be a broad cy- 
clotron feature overlapping the end of the spectrum. 
The fitted blackbody component, however, with kT^B — 
0.641q q4 keV, is typical for an AXP. It has a radius of 
0.8 ds km, and a bolometric luminosity of 1.4 x 10 34 d\ 
ergs s -1 . 

The flux from the second Chandra observation of this 
source is similar to that of the first; however, we cannot 
use it to address the peculiarities of the spectrum because 
of the greatly enhanced background in CC-mode. 

2.2. HESS Jl 731 - 34 7/G353. 6-0.7 

HESS J1731— 347 is an extended source that coincides 
with a patch of diffuse X-ray em i ssion i n t he ROSAT 
All- Sk y Survey (jAharonian et all l2008bO . iTian et all 
(2008) discovered that the X-ray feature is associated 
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with multiwavelength radio emission that comprises a 
new SNR, G353.6— 0.7. A two-dimensional Gaussian fit 
to the TeV source gives semi-major and semi-minor axes 
of oi = 0.°18 ± 0.°07 and a 2 = 0.°11 ± 0.°03 for the main 
part of the e mission, after allowing f or the instrumen- 
tal resolution (jAharonian et al.1l2008b[ ). The TeV source 
thus fits within the 0.°5 diameter radio shell (Figure [Q), 
while the diffuse ROSAT emission is a half shell match- 
ing the southern Galactic half of the radio remnant (see 
ITian et al.l 120081 ). The other half is likely absorbed in 
soft X-rays by the increasing molecular gas density to- 
w ard the Gal a ctic p lane . 

ITian et all (|2008f) assumed that the SNR is associ- 
ated with a n an adjacent H II region at d = 3.2 kpc. 
ITian et ail (|2009f ) found a molecular cloud at the same 
velocity as an H I absorption feature that they con- 
sider to be associated with the H II region. Consid- 
ering the implied ~ 28 pc diameter of the SNR, they 
suggested that G353.6— 0.7 is ~ 27,000 years old and 
has entered the radiative phase. This, they argued, 
is the best evidence that an old SNR can still accel- 
erate part icles to greater t han T eV energies, as pre- 
dicted by lYamazaki et al.l (|2006f ) . and that the TeV 
emission comes from the interaction of the SNR shock 
with the molecular cloud. However, while the outline of 
HESS J1731— 347 follows the radio and X-ray contours, 
it is not entirely clear whether the the TeV emission is 
peaked in the center or follows the shell. So me evidence 
for the latter is seen in newer HESS data (fA cero et &U 
120091 ). There is also a "tail" of TeV emission to the 
west, beyond the radio shell, which may or may no t 
be part of the same source (|Aharonian et al.1 12008b) . 
The energy flux of HESS J1731-347 above 0.5 TeV is 
w 3.3 x 10 -11 erg cm~ 2 s _1 , corresponding to luminosity 
L 7 = 4.0 x 10 34 dl 2 erg s -1 . 

Suzaku, XMM- Newton, Chandra, and Swift observa- 
tions of HESS J 173 1-347 are listed in Table [Q Figured] 
shows the XMM-Newton image, which reveals the fila- 
ments of G353.6— 0.7, as well as a point source near the 
center of the radio shell that was not detected in the 
ROSAT All-Sky Survey. We began by fitting the spec- 
trum of a region of bright, diffuse SNR emission in the 
Suzaku data. Photons from an elliptical region « 6' long, 
adjacent to the point source and extending to the north- 
east of it, were extracted from all three XIS detectors and 
fitted to an absorbed power-law spectrum. The result- 
ing photon index is T = 2.54 ± 0.15, with an absorbing 
column density Ah = (1.65 ± 0.17) x 10 22 cm" 2 . There 
is no evidence for emission lines in the spectral residuals. 
Mor e complete analys i s of the spatially resolved spec- 
tra (|Acero et al.l 120091: ITian et al.l 120091 ) reveal ubiqui- 
tous non-thermal emission, with a gradient of Ah across 
the remnant that is consistent with additional interven- 
ing material toward the Galactic plane. The spectra thus 
indicates the essentially nonthermal nature of the SNR, 
similar to the brightest TeV SNR RX J1713. 7-3946. 

We next turned to the compact source 
XMMU J173203. 3-344518. It lies off-axis in both 
the XMM-Newton and Chandra images. The best posi- 
tion, from the Chandra image, is R.A. = 17 h 32 m 03. s 41, 
deel = -34°45'16."6 (J2000.0). It lacks an optical 
or IR counterpart in the Digitized Sky Survey or in 
2MASS, and is therefore a candidate for a neutron star 



TABLE 3 

Spectral Fits to XMMU J173203. 3-344518 



Parameter 


Suzaku 


XMM-Newton 


Power-law Model 


Nn (10' Z2 cm" 2 ) 

r 

^(0.5 - 10 keV) a 
-M0.5-10keV) b 


3.3 ±0.2 
4.5 ±0.2 
3.5 x lO" 12 
1.8 x 10 35 dl 2 
0.98(136) ' 


3.2 ±0.1 
4.7 ±0.1 
2.8 x 10" 12 
2.0 x 10 35 dl 2 
1.87(148) ' 


Blackbody Model 


jV H (1022 cm" 2) 
kT (keV) 
R (km) 

Fx (0.5 - 10 keV) a 
L x (0.5- 10 keV) b 
Xi(") 


1.52 ±0.09 
0.50 ±0.01 

l.lrfs.2 
3.2 x 10~ 12 
9.1 x 10 33 dl 2 
1.0(136) 


1.40 ±0.06 
0.49 ±0.01 

1.0d 3 .2 
2.8 x 10~ 12 
7.7 x 10 33 dl 2 
1.15(148) ' 


Two Blackbody Model 
,r / 1r ,22 -2> -i on+0.27 -, cc +U.2y 



fcTi (keV) 
iJi (km) 
kT 2 (keV) 
R 2 (km) 

^(0.5 - 10 keV) a 
-M0.5- 10 keV) b 



0.42 



0.18 

+0.07 
0.05 
1.6d 3 . 2 
S5+ - 25 

0.17d 3 .2 
3.4 x 10" 12 

1.1 x 10 34 d§ 2 
0.81(134) ' 



0.40Z 



-0.12 
-,+0.03 
-0.08 
1.5d 3 . 2 

0.64 ±0.10 

0.36 d 3 . 2 
2.8 x 10" 12 
9.4 x 10 33 dl 2 
0.99(145) ' 



a Absorbed flux in units of erg cm -2 s — 1 , average of all detec- 
tors. 

b Unabsorbed luminosity in units of erg s — 1 . 



in G353.6-0.7. 

Power-law spectral fits to the Suzaku and XMM- 
Newton data for XMMU J173203.3-344518 have very 
steep slopes, r = 4.5 — 5.7 as listed in Table [3] and the 
power-law fit to the XMM-Newton data is poor. The 
Chandra spectrum is not as steep (r = 3.98 ± 0.13), but 
it should be disregarded. As a strong source with count 
rate 0.21 s _1 , it is affected by pileup, which tends to 
flatten the spectrum. A better and more physically mo- 
tivated fit to the Suzaku and XMM-Newton data is ob- 
tained with a blackbody of kT « 0.5 and radius of 1 km. 
This is typical for a CCO, but also for an AXP. A good 
fit is also obtained using two blackbodies, with kT± 
0.4 keV and kT 2 in the range 0.6-0.9 keV. The single and 
double blackbody fits are shown in Figured] The higher 
temperature component is also seen in AXPs, but is 
somewhat higher than has been s een in any CCO, which 
tend to cluster around 0.4 keV (IKargaltsev et al.l I2002t 
Lazendic et al.ll2003Hde Luca et aLir2 004; Ha lpern et al.l 
20071 iGotthelf fc Halpernl 120091 : iPavlov fc Lunal 12009) . 



Adopting the two-blackbody model, we obtain a bolo- 
metric luminosity of ~ 1 x 10 34 dl 2 er g s_1 j which is in 
the rang e of quiescent AXPs, but on the high side for 
a CCO (|Gotthelf fc Halpernl 12001 . The fitted column 
densities in the thermal models of the compact source 
are consistent with that of the diffuse nonthermal SNR 
emission, which favors these models and lends further 
support to the association of the compact source with 
the SNR. 

The Chandra image shows faint, symmetric nebulos- 
ity around XMMU J173203.3-344518. Excess emission 
compared to the point spread function, modelled at the 
off-axis location of the piled-up point source, is seen out 
to ~ 50" (Figure [5]). It is not immediately obvious 
whether this is a PWN or a dust scattered halo, although 
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Fig. 4.— (a) Spectrum of XMMU J173203.3-344518 from the 
XMM-Newton observation, fitted to a blackbody model. The upper 
spectrum is from the pn CCD, and the lower spectrum is from the 
two MOS detectors, (b) The same spectra fitted to a two-blackbody 
model (see Section 2.2). Fitted parameters are listed in Table|3] 

its symmetry suggests the latter. We extracted counts 
from an annulus 5.5" < r < 50", and fitted its spectrum 
to a power-law model with column density fixed at the 
same value as for a power-law fit to the point source. 
This yields an acceptable fit with T = 4.97 ± 0.14 and 
a flux that is only 9% of that of the point source. The 
measured spectrum is anomalously soft for a PWN, but 
is slightly softer than that of the point source, which is 
expected for a dust scattered halo. To further check this 
hypothesis, we fitted the spectrum from a smaller annu- 
lus, 5.5" < r < 20" containing the brightest part of the 
halo, and find a spectrum that is consistent with that of 
the point source. In any case, the absence of a brighter 
PWN probably limits the spin-do wn power of th e puta- 
tive pulsar to < 4 x 10 36 erg s" 1 (|Gotthelfll2002T ). 

At the current stage of the analysis there is no strong 
evidence for variability of the point source. The Suzaku, 
XMM-Newton, and Chandra observations agree to within 
20% in flux. Two short observations with the Swift 
XRT listed in Table Q] and separated by 39 days are 
useful to search for X-ray variability that would be 
a signature of a magnetar. The 0.5-10 keV counts 
from XMMU J173203. 3-344518 are marginally consis- 
tent between these two observations, with count rates 
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Fig. 5.— Nebulosity associated with XMMU J173203.3-344518 
in the Chandra image. The radial profile (histogram) is compared 
to a CHART/MARX simulation of a point source (solid curve) using 
the position and spectrum of XMMU J173203. 3-344518, and al- 
lowing for CCD pileup. "Field star" refers to the point source to 
the southeast. Inset: The Chandra image in which the compact 
component of XMMU J173203. 3-344518 has been excised (solid 
circle). The dahed circles, of radii 20" and 50", enclose the regions 
used for the spectral fits and flux calculations discussed in Section 
2.2. 

of 0.056 ± 0.006 s- 1 and 0.073 ± 0.007 s" 1 , respectively. 
The significance of the difference is < 2 a. These rates 
are consistent with the fluxes measured previously for the 
source by Suzaku and XMM-Newton. Therefore, no evi- 
dence of magnetar-like activity has been seen over the 2 
year time span during which XMMU J173203.3-344518 
has been observed. 

During the XMM-Newton observation, the EPIC pn 
CCD was operated in full frame mode with a time res- 
olution of 73.4 ms. This allows us to search for periods 
> 150 ms. After barycentric correction, we find only 
marginal evidence for a period of 1.010441(6) s from 
XMMU J173203.3-344518, modulated with a - 10% 
pulsed fraction. For this signal, Z\ = 32.7, correspond- 
ing to a 0.4% chance probability for a blind search for 
P > 1 s, or 2.7% chance probability for a search of all 
independent periods > 150 ms. We consider this too 
weak to claim a secure detection, although such a period 
is plausible in either the AXP or the CCO scenario, the 
implications of which are discussed below. 

3. DISCUSSION 

Previous authors have argued that particles accel- 
erated at the SNR shock are responsible for the TeV 
emission via pion decay from both HESS J1713— 381 
dAharonian et alj l2008aft and HESS J1731-347 
( Tian et alj I2008T ). In the latter case, interaction 
with nearby molecular clouds is envisioned (jTian et alJ 
120081 I2009t i. On the other hand, neither TeV source 
clearly has a pure shell-like morphology, so it cannot 
be ruled out that their associated pulsars make some 
contribution to the TeV flux via the inverse Compton 
mechanism. 

Magnetars are not expected to power TeV PWNe be- 
cause their large magnetic fields have already spun them 
down to where their present spin-down power E is in 
the range 10 32 — 10 35 erg s^ 1 . This is much less than the 
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E = 10 36-38 erg s _1 of pulsars powering TeV sources and 
PWNe in general. Also, unlike ordinary pulsars, it is not 
certain that magnetars can accelerate particles to TeV 
energies. Magnetar models that involve strong currents 
on closed, twiste d magnetic field lines develop vo ltages 
of only - 10 9 V (|Beloborodov fc Thompson] 12001 . On 
the other hand, it is not excluded that the ordinary pul- 
sar mechanism operates on open magnetic field lines of 
magnetars, which may develop particle-dominated winds 
that flow out to become shocked PWNe. Evidence for 
this possibility comes from the two magnetars that are 
also transient radio pulsars ( Camilo ct al. 2006, 2007). 

Even though no TeV source has yet been associated 
with a known magnetar, we consider it possible that 
a magnetar could energize a PWN at a younger stage, 
when its period is short and its spin-down power is 
large. Then, some high-energy electrons may escape 
the high .B-field region and emit via inverse Compton 
scattering into the TeV band for an extended time. A 
possible prototype of such systems is PS R J1846— 0258, 
the " transitional" ~ 700 yr old pulsar (Gotthelf ct al. 
l2000h associat ed with HESS J1846 -029 in the shell-like 
SNR Kes 75 ()Diannati-AtaJl2008fh This 0.32 s pulsar 
has a dipole field that approaches magnetar strength 
(B s = 4.9 x 10 13 G) and displays AXP-like bursts (Gavriil 
et al. 2008) ; it is therefore likely of an intermediate class 
connecting the rotation-powered and the magnetar pul- 
sars. It is also one of the most energetic pulsars; its 
E = 8.1 x 10 36 erg s _1 is easily sufficient to power its 
X-ray PWN and the TeV emission of HESS J1846-029. 

The quantitative details of any magnetar source model 
for extended TeV emission can be constrained by inde- 
pendent age estimates from associated SNRs. The num- 
ber of magnetars securely associated with shell SNRs 
is arguably at most three, IE 2259+586 in CTB 109, 
IE 1841-045 in Kes 73, and SGR 0526-66 in the LMC 
SNR N49 (Gaensler et al. 2001). To these we may 
tentatively add CXOU J171405. 7-381031 in CTB 37B, 
and IE 1547.0 —5408 because of its possi ble radio shell 
G327.24-0.13 (|Gelfand fc Gaenslen 12001 . In addition 
to providing age estimates, such associations are valu- 
able for inferring properties such as natal kick velocity, 
and energetics of a supernova that is able to create a 
magnetar (e.g., Vink & Kuiper 2006). The radio emit- 
ting AXP XTE J1810-197 is the only magnetar that 
has a well measured pro per motion, corresp onding to 
v t = 212 ±35 d 3 . 5 km s" 1 (jHelfand et al.|[2007h . which is 
not unusual compared t o ordinary young neutron stars. 
lAharonian et a l. (2008a) have already pointed out that 
the displacement of CXOU J171405. 7-381031 from the 
center of the shell of CTB 37B would require a trans- 
verse velocity of ~ 1000 km s _1 to the east for an age 
of 5000 yr. The velocity would be even larger if an as- 
sociation with th e historical supernova of 373 AD is as- 
sumed. However. lAharonian et all (|2b08a) also note that 
the brighter X-ray and radio emission on the eastern side 
(see Figure [lj could indicate a higher density density 
there, implying that the pulsar is closer to the explosion 
center, reducing its inferred velocity. 

In the case of XMMU J173203.3-344518, the lack 
of secure spin measurements means that the alterna- 
tive of a low magnetic field neutron star is not yet 
ruled out. Its tentative 1 s period is intermediate be- 



tween those of the magnetars (2-12 s), and the cen- 
tral compact objects (CCOs) in SNRs that manifestly 
have weak fields and spin periods between 0.1 and 0.4 s 
(|Gotthelf fc Halpernl[200l . CCOs are not clearly distin- 
guished from AXPs on the basis of their X-ray spectra 
alone, in the absence of spin-down data. The temper- 
ature^) and X-ray luminosity required to fit the spec- 
trum of XMMU J173203. 3-344518 are similar to those 
of CCOs, which require either one or two blackbodics 
of kT between 0.2 and 0.6 keV, and 10 33 < L x < 10 34 
erg s . However, a fraction of magnetars are observed 
for many years in quiescent states of 10 33 < L x < 10 34 , 
with spectra that are similar to CC Os. Examples are the 
trans ient AXPs XTE J 1810-197 (I Gotthelf fc Halpernl 
|200"1 . IE 1547.0-5408 dGelfand fc Gaenslerl 120071) . and 
CXO J16 4710.2-455216 (iMuno et all 120061 and SGR 
1636-41 (jMereghetti et al.l 120061 Because such low- 
luminosity sources are difficult to classify before they 
are seen in outburst, they may represent the majority of 
magnetars. Eventually, magnetars show X-ray variabil- 
ity, whereas weak jB-field CCO pulsars have not. The 
classification of XMMU J173203. 3-344518 remains un- 
certain until either X-ray variability is established or its 
spin parameters are measured. 

4. CONCLUSIONS 

We have identified two candidate magnetars in SNRs 
based on their X-ray spectra, probable variability, and 
in one case, CXOU J171405. 7-381031 in CTB 37B, 
a typical magnetar pulse period of P = 3.82 s. 
We have a second timing observation of each candi- 
date planned that will measure its all-important period 
derivative, from which the spin-down luminosity, char- 
acteristic age, and dipole magnetic field strength will 
be derived. This will test the magnetar hypothesis for 
CXOU J171405. 7-381031. Evidence of a hard spectral 
excess also needs to be confirmed and explored at higher 
energy. 

In the case of XMMU J173203. 3-344518 in 
G353.6— 0.7, only marginal evidence of a 1.01 s pe- 
riod is found. However, given its soft X-ray spectrum 
typical of magnetars, a period such as this is expected 
and is of great interest. It lies between the low end 
for magnetars (IE 1547.0-5408; P = 2 s) and the 
transitional object PSR J1846-0258 (P = 0.32 s). 
With such a short period, it is possible that the E of 
XMMU J173203.3-344518 exceeds the TeV luminosity 
of HESS J1731-347, « 4.0 x 10 34 erg cm" 2 s"\ or did 
so in the recent past. Even if this candidate rotation 
period is not confirmed, it is possible that the true 
period will be discovered in more sensitive observations 
of XMMU J173203.3-344518 that will establish it 
as a magnetar. On the other hand, if further spin 
measurements show that XMMU J173203.3-344518 
is a low B- field CCO, then its spin-down power was 
never sufficient to power HESS J1731— 347. In that 
case, the SNR G353.6— 0.7 must be responsible for 
the TeV emission. Determining the spin properties of 
XMMU J173203. 3-344518 is therefore a crucial step in 
testing whether old SNRs can power TeV sources. 
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